Vol. 97 (2000) ACTA PHYSICA POLONICA A No. 1(2) 



Proceedings of the European Conference "Physics of Magnetism 08" , Poznah 2008 

Zero temperature phase diagram of a d-wave 
superconductor with Anderson impurities 

L. S. BORKOWSKI 

Quantum Physics Division, Faculty of Physics, 
A. Mickiewicz University, Umultowska 85, Poznan, Poland 

(Version date: June 10, 2008) 

We study the model of a d-wave superconductor interacting with finite 
concentration of Anderson impurities at zero temperature. The interaction 
between impurity and conduction electrons is taken into account within the 
large-A'^ approximation. We discuss the obtained phase diagram and its 
dependence on the main energy scales. 

PACS numbers: 74.81.-g,74.25.Dw 

Magnetic and nonmagnetic impurities in correlated electron systems are 
an important probe of the properties of the host material. There is large and 
growing body of research on impurities in high-temperature superconductors. [1, 2] 
The refinement of experimental techniques probing local properties stimulated the 
theoretical work on defects in those systems. The investigation of the effects of Zn, 
Ni and other dopants in YBCO and BSCCO answered some important questions 
and raised new ones. 

Measurements on some compounds show that the superconducting order 
parameter is not uniform over the entire sample. [3]- [7] Scanning tunnelling spec- 
troscopy measurement showed that modulation of the structure of Bi2Sr2CaCu208+x 
is correlated locally with the magnitude of the energy gap. [8] The spatial varia- 
tion of Ao(r) may result e.g. from the structural supermodulation affecting the 
strength of local pairing interaction. [9] 

The intrinsic spatial variation of the superconducting gap raises the possi- 
bility of observing impurity states on both sides of the quantum phase transition 
in the same sample. Theoretical work on magnetic impurities in systems with 
reduced density of states near the Fermi surface[ll]-[18] showed that the resonant 
impurity states may be viewed as a sensitive probe of the superconducting state. 
If the coupling to the magnetic impurity is small compared to the energy scale 
associated with the gap, the impurity is decoupled from conduction electrons. 
This impurity quantum phase transition occurs at finite coupling, provided the 
particle-hole symmetry is broken, and may be studied by STM techniques. The 
low-energy behavior of the model depends on the exponent r in the conduction 
electron density of states, N{e) ~ |e|'', where the Fermi level is fixed at e = 0. 
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The interaction of impurities with the conduction electron band may be 
studied in the Anderson model with a BCS-type pairing interaction, 
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This model allows studying also the mixed valence regime where the impu- 
rity occupation number is less than 1 and charge fluctuations are dominant. We 
assume a two-dimensional d-wave order parameter of the form A(A:) = Aq cos(2^), 
where (j) is the angle in the kx — ky plane. The constraint of single occupancy of 
the impurity site by adding a term A(n/ — 1) to the Hamiltonian, where A is the 
Lagrange multiplier and taking A — > oo. Minimizing the free energy with respect 
to the resonant level energy e/ and z =< >=< b > and taking the mean field 
approximation we obtain 
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Equations equations (2) and (3) are solved self-consistently with the gap 
equation 
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The self-energy in the full conduction electron Green's function G is av- 
eraged over impurity positions. We solve these equations self-consistently and 
obtain the phase diagram as a function of the parameters of the model. 

The obtained phase diagram is shown in Fig. 1. The approximate slope 
of the impurity quantum phase transition line, —0.26, agrees with results of the 
numerical renormalization group (NRG) method[15] and large-A?^ one-impurity 
calculation. [18] 

When the bare impurity level Eq lies closer to the Fermi energy the self- 
consistent treatment of finite concentration of impurities leads to reentrant behavior. [19] 
In that part of the phase diagram the pair-breaking is weaker and is spread over 
wide energy range. 

Fig. 2 shows the impurity occupation number at the impurity decoupling 
transition. It reaches 1 at the transition. For Fq 3> Ao, nj(i?o) slowly approaches 
1, |<iri//(i£^olro=const ^ as E^o/Fq (-^o/ro)CTitical- However for small Fq the 
dependence of n/ on Eq/Fq becomes singular and the slope |<^?^//c^£'o|ro=const ^ 
1. 
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Fig. 1. The phase diagram of the d-wave superconductor with nondegerate Anderson 
impurities for several impurity concentrations. The lines are guide to the eye. The 
order parameter amplitude is Ao = O.OID, where D is half of the conduction electron 
bandwidth. All energies arc scaled in units of D. The dotted line indicates the impurity 
quantum phase transition. 
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Fig. 2. The impurity occupation number as a function of the ratio Eq/Fq for several 
values of Fq. At the impurity transition Uf 1. 
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Fig. 3. The location of the critical point (dotted line) and the impurity decoupling 
transition (broken line). 

In the limit of vanishing interaction Fq — > 0, E'o = is the singular point of 
the model. For £^o > the superconducting state survives even for large impurity 
concentration. 

For any finite n there is a critical point (Foc-Eoc)- For Fq slightly larger 
than Foc there is a superconductor-normal state transition at Eqi < Eqc and 
another normal-superconductor transition at Eq2, where £^01 < Eq2 < Eqc- At 
fixed n, Eq2 — Eqi ~ a(n)(Fo — Fqc), where a(n) weakly depends on n. The 
location of this critical point in the Eq-Tq plane is shown in Fig. 3. 

It would be interesting to test this theoretical picture in experiment. Near 
the impurity quantum phase transition the impurity state is very sensitive to 
small changes of hybridization Fq or impurity level energy Eq. In compounds 
with spatially varying energy gap this could lead to impurities existing on the two 
sides of the transition line in different parts of the sample. 

The phase diagram calculated in this work might also be relevant in some 
heavy-fermion compounds where similar competition occurs between energy scales 
associated with the Kondo screening and the superconducting correlations. Stud- 
ies of CeCu2(Sii_xGex)2 under varying hydrostatic pressure reveal two supercon- 
ducting domes in the phase diagram. [20, 21] The existing interpretation of this 
dependence on pressure relies on additional valence-fiuctuation mediated pairing 
mechanism. [22] 

However our work suggests that the second stiperconducting dome in 
CeCu2(Sii_xGex)2 at high pressure may follow from weakened pair-breaking. The 
change of pressure shifts the chemical potential and brings the system to the 
mixed- valence regime when the bare /-level E^ of Ce ions approaches Ep- The 
phase diagram in Fig. 1 shows that in this limit superconducting correlations are 
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less affected by pair-breaking. 

The large- method used in the present calculation gives qualitatively sim- 
ilar results for larger N. Also the symmetry of the order parameter should not 
introduce drastic changes to the phase diagram. The reentrant behavior as a func- 
tion of Eq results from the competition between the formation of the impurity 
resonance and superconducting correlations and depends mainly on the ratio of 
the relevant energy scales. 

Extension of the theory beyond the mean field should not change the phase 
diagram qualitatively. A more detailed description of physics in the vicinity of 
the impurity transition line requires careful treatment of low-energy scattering in 
specific superconducting compounds. 
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